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ABSTRACT 

Active Galactic Nuclei (AGN) show excess X-r ay em ission above 10 keV compared with ex- 



trapolation of spectra from lower energies. Risaliti et al.l have recently attempted to model the 
hard X-ray excess in the type 1.8 AGN NGC 1365, concluding that the hard excess most likely 
arises from Compton-scattered reflection of X-rays from an inner accretion disk close to the black 
hole. Their analysis disfavored a model in which the hard excess arises from a high column density 
of circumnuclear gas partially covering a primary X-ray source, despite such components being 
required in the NGC 1365 data below lOkeV. Using a Monte Carlo radiative transfer approach, 
we demonstrate that this conclusion is invalidated by (i) use of slab absorption models, which 
have unrealistic transmission spectra for partial covering gas, (ii) neglect of the effect of Comp- 
ton scattering on transmitted spectra and (iii) inadequate modeling of the expected spectrum of 
scattered X-rays. The scattered spectrum is geometry dependent and, for high global covering 
factors, may dominate above lOkeV. We further show that, in models of circumnuclear gas, the 
suppression of the observed hard X-ray flux by reprocessing may be no larger than required by the 
'light bending' model invoked for inner disk reflection, and the expected emission line strengths 
lie within the observed range. We conclude that the time-invariant 'red wing' in X-ray spectra is 
probably caused by continuum transmitted through and scattered from circumnuclear gas, not 
by highly redshifted line emission, and that measurement of black hole spin is not possible. 
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X-ray observations of AGN reveal evidence for 
absorbing g as surrounding an accre ting black hole 
system (see Turner &: Miller 20091 for a review) 
and UV spectroscopy si milarly shows comple x 



absorption systems (e.g. ICrenshaw et al.l 120031 ) 



Fe XXV and Fe xxvi absorption lines provided 
the first clear evidence for gas with hydrogen 
column densities A'h ^ ^0^'^^^^ cm^^ and out- 
flow velocities up to 0.35c may be associated with 
highly-ionized zone s in up to 40% o f type 1 AGN 
(|Cappi et all l2009t iTombesi et all l2010l 12011 



20121 ). A recent Suzaku study (jTatum et al.ll2013D 
of a hard X-ray selected sample of type 1 AGN 
has shown luminosities above lOkeV to be much 
higher than expected from extrapolation of mod- 
els from the 0.5-10 keV band. This 'hard excess' is 
ubiquitous in local type 1 AGN and, owing to the 
sharp absorption edges associated with the hard- 
est source spectra, cannot be explained other than 
by the presence of a Compton-thick layer of low- 
ionization absorbing gas, covering a large fra ction 
of the continuum source ( Tatum et al.ll2013 ). 
Partial-covering absorption from X-ray absorb- 
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ing gas below the Compton-thick limi t has long 
been known in local type 1 AGN (jPiro et al 



NuSTAR, and for Compton -thick gas, it i s cru - 



20051 : iMiller et alJ 120081 ) and changes in the ab- 



sorber provide a natural explanatio n for spectral 
variab ility on long timescales. The iTatum et al 



( 2013f ) analysis extends the parameter space for 
partial-covering absorption into the Compton- 
thick regime. 

It has long been claimed that X-ray spectra also 
show features caused by Compton scattering and 
absorption (known as 'reflection') in the outer lay- 
ers of an inner accretion disk, blurred by relativis- 
tic effects. Spectral curvature over the 2-8 keV 
band may be fit with a model that convolves re- 
flection spectra (.Ross fc Fabian .2005i) with gen- 
eral relativistic effects fc.f. lLaoiill99l[ ). 

However, it is difficult to distinguish absorption 
and reflection signatures based solely on fitting 
individual X-ray spectra, because Compton scat- 
tering and absorption shape both models, leading 
to similarities in overall shape. Consideration of 
AGN variability characteristics can help, and per- 
haps the biggest problem for the blurred reflection- 
dominated picture is the lack of any clear corre- 
lation between continuum and reflection flux in 
variable AGN (sectiong]). 

A particularly interesting case is the type 1.8 
AGN NGC 1365, which has b een observed exten- 
sively in the . 5-10 ke V range ()Risaliti et al.ll2000l 
2005bllal. l2007i |2009cP and at higher energies by 



Risaliti et al.l ( 2009a , 2013 ). These observations 



have shown spectral variability on timescales as 
short as a few hours, attributed to changes in the 
X-ray absorber. 



Risaliti et al. ( 2013t ) present simultaneous XMM- 



Newton and NuSTAR data for NGC 1365, cover- 
ing 0.5-79 keV. The authors claim that the data 
show relativistic disk features: broad Fe Ka line 
emissio n and Compto n scattering excess above 



lOkeV. iRisaliti et al.l claim that the reflected 
component arises within 2.5 gravitational radii of 
a rapidly spinning black hole and that absorption- 
dominated models that do not include relativistic 
disk reflection can be ruled out. As this is a very 
strong assertion, we briefly review the published 
model upon which it is based (Risaliti et al. 2013). 

There are two key processes that shape AGN 
X-ray spectra: Compton scattering and photoelec- 
tric absorption. At the high energies observed by 



cial to inclu de both effects. Risaliti et all (l2013h 



used XSTAR ([Kallman fc Baut ista 2001) to create 
a model representing a slab of ionized absorbing 
gas. However, the fits for NGC 1365 required col- 
umn density iVn ~ 5 x lO^^'cm^^. Such gas would 
be Compton-thick, but xstar does not include 
Compton scattering and is not a valid model for 
column densities in this regime (in addition, we 
note that currently-released versions of xstar do 
not calculate opacity above 20 keV) . 

Even if Compton scattering were included, the 
effect on t he spectrum d e pends strongly on gas 



geometry. IRi saliti et al. (2013) considered two 



'bounding possibilities' for gas geometry, but we 
shall show below that these estimates are not ap- 
propriate. Finally, the authors used another ab- 
sorption model, PLCABS, to estimate the Comp- 
ton scattering correction they should apply to the 
observed luminosity. Not only is such a simple 
flux- loss consideration inadequate (see section[3]), 
we note that plcabs is only valid up to lO keV for 
such high column densities ( YaaoobI 19971 ). 

In this paper we present Monte Carlo simu- 
lations that demonstrate the importance of the 
inclusion of Compton scattering for the column 
density regime apphcable to NGC 1365. We also 
demonstrate the strong geometry-dependence of 
the predictions and discuss the implications for 
modeling and interpretation of AGN X-ray spec- 
tra. 

2. Absorption and scattering by a spheri- 
cal cloud 

2.1. The importance of geometry 

We first consider absorption and scattering by 
a single cloud of constant density, absorbing pho- 
tons from a distant source, to demonstrate the 
key effects that must be taken into account. The 
as is assumed to be cold, with solar abundances 
1989| ) and cross-sections of 



Anders fc Grevesse 



Verner et al.l ( 1996[ ). We first consider the effect 



of cloud geometry alone, by omitting the effects 
of Compton scattering, for a cloud with column 
densit y 5 x lO^^'cm"^, as found by IRisaliti et al. 
(|2013l ). The cloud is illuminated by a power- 
law of photon index F = 2 and has a normaliza- 
tion, that would be seen by the observer in the 
absence of absorption, of unity at 1 keV. Fig.[T] 
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plots the transmitted spectrum Ef{E), for pho- 
ton energy E and energy flux per unit energy 
interval f{E). The lowest dashed curve shows 
the standard calculation obtained for a plane slab 
of uniform density, where the transmitted flux 
f{E) oc exp{—Niia{E)) for photoelectric absorp- 
tion cross-section per hydrogen atom (j{E). 

However, a plane slab is an unlikely choice for 
a clumpy distribution of gas that partially covers 
the source. The distribution of column density is 
likely to be complex, and it is instructive to see 
the effect of changing the geometry to that of a 
uniform spherical cloud. Now, there is a distribu- 
tion of column density associated with the cloud, 
p(7Vh) = 27VH/iV2_^^, for Nn < iVn^max, where 
-^H,max is the maximum column density through 
the cloud. Integrating over this distribution, and 
in the absence of Compton scattering, we expect 
transmitted flux 



2 (1 - [1 + iVH,maxg(^)] e- 



(1) 

The upper dashed line in Fig. [1] shows the results 
for a cloud with mean density the same as the 
slab and for which A^H.max is 3/2 times the mean 
column density. At high optical depths the trans- 
mitted flux f{E) — > 2 (iVH.maxfl-E))"^, and since 
a{E) has an approximately power-law form, a 
hard, powerlaw spectrum develops. In this regime, 
the dominant contribution to the transmitted flux 
comes from the periphery of the cloud, rather like 
the 'silver lining' effect seen in illuminated terres- 
trial clouds. 

The spherical cloud geometry has a significantly 
different transmitted spectrum from the slab ge- 
ometry. Such a cloud model has the same num- 
ber of model parameters as the slab model, and 
is likely to be a better representation of a partial- 
covering cloud than a plane slab. This calculation 
alone demonstrates how careful one must be when 
fitting absorption models to AGN X-ray spectra. 

2.2. The effect of Compton scattering 

At high column densities, incident radiation is 
Compton scattered out of the line of sight. To 
model the transmitted and scattered spectra, we 
carried out a Mon te Carlo radi a tive t ransfer sim- 
ulation, following iMiller et all (|2009l ). Muhiple 



absorbed by a K shell of Fe, an Fe K l ine pho- 
ton could be produced ( Miller et al. 20091 ). Angle- 
dependent Klein-Nishina Compton scattering was 
calculated. The Monte Carlo spectra have been 
smoothed with a gaussian of FWHM AE/E = 
0.0283, appropriate for the resolution of typical 
CCD spectral data around FeK. 

The transmitted spectrum is shown by the 
lower, solid, red line in Fig.[T] At high energies 
the transmitted flux is significantly attenuated by 
Compton scattering. However, the spectral shape 
is significantly different from either the slab or 
sphere when Compton scattering is neglected. At 
low energies, the transmitted flux is again dom- 
inated by light leaking through the periphery of 
the cloud, so that the spectrum tends towards the 
same form as the spherical, non-scattering cloud. 

However, the light that is lost due to Compton 
scattering reappears in the scattered spectrum, 
shown in Fig.[T]by the upper solid curve. The spec- 
trum has been integrated over all solid angles, so 
that the normalization is such that this is the scat- 
tered spectrum that would be measured by an ob- 
server if the primary source were fully surrounded 
by such clouds, but neglecting the effect of radia- 
tion being reabsorbed or scattered by other clouds 
(which will be discussed in section[3]). If instead 
only some fraction of lines of sight are intercepted 
by clouds, the normalization should be multiplied 
by the global covering factor. For any global cover- 
ing factor Cg ^0.1, the scattered spectrum dom- 
inates over the transmitted spectrum, and both 
spectra have considerably different shapes from 
the non-scattering slab and sphere. It is worth 
noting th e low equivalent width of the Fe Kg ern is- 
sion line (|Miller et al.ll2009l: lYaqoob et al.ll2010t ). 

Finally, Fig.[T] shows a comparison with reflec- 
tion from a plane, neutral, optic ally-thick slab as 

calculated by the PE XRAV model ( Magdziarz fc Zdziarskil 
19951 ) within XSPEC ()Arnaudlll996l ). without calcu- 



scatterings were allowed, and when photons were 



lation of the emission line strength. The normal- 
ization and shape of the PEXRAV spectrum varies 
within inclination: the curve shows the case for 
cos 9 = 0.5, which matches most closely the cloud 
scattered spectrum around lOkeV. The PEXRAV 
spectrum is qualitatively similar to the cloud spec- 
trum, although there are differences in shape, es- 
pecially above 20keV when the cloud spectrum 
has a less pronounced 'Compton hump'. Given 
also the significant differences in the transmitted 
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Fig. 1. — Transmitted and scattered spectra 
from a single cloud: lower, dashed: transmitted 
spectrum through a uniform plane slab, neglect- 
ing Compton scattering; upper, dashed: transmit- 
ted spectrum through a uniform spherical cloud 
of same mean density as the plane slab, neglect- 
ing Compton scattering; lower, solid: transmit- 
ted spectrum through a uniform spherical cloud, 
including Compton scattering; upper, solid: scat- 
tered spectrum from the uniform spherical cloud; 
dotted: reflection from a plane slab calculated by 

PEXRAV. 



spectra, we cannot expect to model even simple 
Compton-thick clouds with combinations of ab- 
sorption by plane slabs and planar reflection. We 
shall see in the next section that, in models with a 
high covering factor of clumpy circumnuclear gas, 
the scattered spectrum shows much larger depar- 
tures from the PEXRAV model. The importance 
of geometry for the observed sp ectrum has pre- 
viousl y been demonstrated by M urphv &: YaqoobI 



br a smooth toroidal reprocessor. 



3. Scattering and absorption by a complex 
distribution of circumnuclear gas 

In this section we consider the case of a high 
global covering factor of clumpy gas, when scatter- 
i ng of photons between c louds cann ot be ignored 



As 



Tatum et al 



([Nandra fc Georgel Il994l ) 
20131 ). a simple, clumpy distribution was cre- 



ated by placing 1000 spherical clouds of equal 
radius between inner and outer radii r^in and 
''max from a central X-ray source. Gas was 
placed inside each cloud with a uniform density 
(the density was not increased within overlapping 
clouds). For illustration, we chose a single value 
of A'h = 1-5 X lO^'' cm^^ per unit cloud-radius, 
chosen to approximately maximize the hard ex- 
cess. In order to simulate the effect that in type 2 
AGN we might be looking through a denser distri- 
bution of clouds than in type 1 AGN, the number 
density of clouds increased towards the equatorial 
plane as sin (/), where is the polar angle of ob- 
servation. In Fig. [2] we show two models: model A 
has rmin = 10, Tniax = 20 in units of the cloud 
radius; model B has rmin = 10, Tmax = 30. The 
fraction of sightlines that were empty of gas was 
0.35 at ~ 0, decreasing to < 0.15 for </> > 60° in 
model A, and 0.55 at ~ decreasing to < 0.25 
for (p > 60° in model B. We emphasize that these 
are simplified models of circumnuclear material, 
designed to illustrate the radiative transfer effects, 
and are not to be taken as a proposal for the actual 
gas distribution. 

Spectra were created using 6 x 10^" input 
quanta with energy 0.1 — 400 keV. The input spec- 
trum was a powerlaw of ph oton index F = 2.2, a s 
inferred for NGC 1365 (e.g. iRisaliti et al.ll2009a[ ). 
It was assumed that photons that propagated 
through the equatorial plane were absorbed by 
an accretion disk and, to avoid possible confu- 
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Fig. 2. — Spectra from the models of circmTmuclear gas in sectionOfor models A (upper panels) and B (lower 
panels). Each panel presents spectra from a different azimuthal angle, selected to show, from left-to-right, 
a low, intermediate and high flux state. The panels show the incident continuum (dashed powerlaw line), 
the spectrum of photons that have passed through holes in the gas distribution (green powerlaw line, not 
present in the left panels), the absorbed spectrum transmitted through the gas (blue curve), the spectrum 
of scattered light, including FeK line emission (red curve), and the total spectrum (black solid curve). The 
left panels also show a pexrav model (dotted curve). 
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sion with disk reflection, the disk did not re-emit 
photons. As in section[2] the gas was assumed 
cold and angle-dependent Klein-Nishina Compton 
scattering was calculated. Spectra can be accu- 
mulated only by averaging over finite solid angle, 
and here we adopted a solid angle corresponding 
to that subtended from the central source by a 
cloud at Tmax- A point source was placed at the 
center of the cloud distribution, but the averaging 
of photons over solid angle has the effect of mim- 
icking the partial covering of an extended source 
of comparable size to a cloud. 

Fig.[2]shows example spectra obtained for mod- 
els A and B, viewed at polar angle (/) = 60°, for 
three azimuthal angles. The spectra have been 
smoothed as in section[2l but Monte Carlo noise 
is visible at low flux values. The spectra cover a 
range similar to that observed in the AGN popu- 
lation, from 'reflection dominated', through inter- 
mediate states where primary continuum is visi- 
ble, to high states where the primary continuum 
dominates. All states have signiflc ant hard ex- 



cesses above 20keV, as discussed bv lTatum et al 
(|2013l ). Unobscured sightlines have total light at 



high energy that exceeds the input continuum am- 
plitude. The panels show the contributions to the 
observed spectrum from various components: pho- 
tons that pass through the holes in the distribu- 
tion, photons that pass through clouds but are not 
scattered, and scattered photons, including Fe Ka 
and FeK/3 line emission. The incident continuum 
is also shown as a dashed line. The chief results 
are as follows. 

(i) The transmitted, unscattered spectrum (blue 
curves) is variable and does not look like the spec- 
trum expected from transmission through a slab 
model (section[2]). 

(ii) The scattered radiation is different from that 
calculated in either the isolated cloud or pexrav 
models (the upper left panel of Fig. [2] shows a 
PEXRAV model for photon index F = 2.2 and 
cos 9 = 0.5). The hardness of the scattered spec- 
trum increases with global covering factor. 

(iii) Within each model, the scattered radiation is 
largely invariant with viewing angle, because it 
does not depend on whether the primary source 
happens to be covered, or not, along any one sight- 
line. 

(iv) In all but the highest flux states, the scattered 
light is the dominant spectral component above 



20 keV, where the total amplitude is only a factor 
3-4 below the incident continuum, even in this 
regime where the 2 — lOkeV continuum may be 
suppressed by much larger factors, 
( y) The line flux o f Fe Ka is low, as expected 
(jMiller et all 120091 : lYaaoob et all l2010t ) . Equiva- 
lent widths vary from 1500 eV in the reflection- 
dominated state down to 10 eV in the highest 
state. The equivalent width of FeKa found in the 
three Suzaku observations of NGC 1365 is 330- 
520 eV (Tatum et al., in preparation). 

4. Discussion 

4.1. The 'red wing' in AGN spectra 

The results of section[3] show that, on long 
timescales where spectral variability may be dom- 
inated by variations in observer covering fraction, 
there should be a spectrum of largely invariant 
scattered light that creates an apparent 'red wing' 
at energies below the Fe K edg e. Such an invariant 
red w ing is well known (e.g. iFabian fc Vaughan 
2003), and using principal co mponents analysis. 
Miller et al.l (|2007l l2008l . hOldi} showed that there 
is an approximately invariant 'offset' component in 
AGN spectra. The effect leads to a problem with 
the explanation of the red wing as arising from re- 
flection by the inner accretion disk. In that model, 
the disk is close to the primary X-ray source, and 
flux variations in the primary source should be 
closely tracked by th e reflected light, but appar - 



ently they are not. iFabian fc VaughanI (|2003f ). 



Miniutti et al. 1 (l2003h and subsequent papers pro- 
posed a 'light bending' model in which a compact, 
primary source near the black hole moves towards 
or away from the disk, such that bending of light 
causes an apparent variation in flux from the pri- 
mary source, b ut not from the refle cted light. The 
calculations of Miniutti fc Fabia"nl (|2004 ) indicate 
that this model has to be finely tuned to obtain 
invariant reflected flux. 

Such an a posteriori model is not required if the 
red wing arises from a combination of transmitted 
and scattered continua from circumnuclear g 
suggested by iTurner et al. ( 2007 ). and not from 
highly redshifted line emission. A lower global 
covering factor than assumed above would pro- 
duce proportionally less scattered light, but the 
general features would be similar to those shown. 
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4.2. Modeling the spectrum of NGC 1365 

In modeli n g the spectrum of NGC 1365, 



Risaliti et al.l (j2013l) recognize tlrat partial cov- 
ering absorption is required, regardless of whether 
there is also a component of blurr ed reflection. 



In an initial comparison (Fig. 3 of iRisaliti et al 



2013I ). a model of blurred reflection plus one ab- 
sorber component was compared against a model 
of two absorber components. The models were 
only fit to the XMM-Newton data below lOkeV, 
and it was then seen how well each model ex- 
trapolated to higher energy. However, the blurred 
reflection included Compton scattering, but the 
two-absorber model did not. Hence it is no sur- 
prise that the first model extrap olated better 



above lOkeV. In a subsequent test, IRisaliti et al 



(|2013l) fit also to the hard X-ray NuSTAR data 
and claimed that the blurred model is statistically 
better. But the transmitted and scattered light 
components were modeled with non-scattering, 
slab absorbers plus a pexrav component, so such 
a test is invalid. We emphasize that the expected 
scattered light spectrum extends below lOkeV 
and affects the XMM-Newton spectrum as well as 
the NuSTAR spectrum, hence the fits to XMM- 
Newton data alone are incorrect if there is high 
column density gas. Point (ii) above demonstrates 
that existing models such as PEXRAV have an in- 
appropriate geometry for the scattered light and 
cannot account for reprocessing by a complex cir- 
cumnuclear distribution. 

The results also demonstrate that estimates of 
the continuum suppression arising from Compton 
scattering (e.g. jlis_aliti et aL 2013) require scat- 
tered light to be properly included. The actual 
suppression obtained, a factor 3-4 at high en- 
ergy, is comparable to the fraction of light lost in 
the 'light bending' model of inner disk reflection, 
where up t o 75% of photon s may be lost into the 
black hole (iMiniutti fc Fa bian 2004). 



4.3. 



Towards improved models of absorp- 
tion and scattering 



The models presented here have assumed gas to 
be cold, whereas in practice we expect it to be ion- 
ized. The models illustrate the importance of the 
effects of geometry and Compton scattering, but 
cannot be used to make detailed fits to AGN spec- 
tra. Ionization balance is calculated by xstar. 



but assuming spherical symmetry and neglecting 
Compton scattering. A promising approach to 
calculating the full 3D radiative transfer through 
ionized circumnuclear gas has been r nade using a 
Monte Carlo approach bv lSim et all (|2008l[201o[) . 
but those models need to be generalized to clumpy 
gas distributions. 

5. Conclusions 

AGN spectra may contain significant levels of 
transmitted and scattered continua from partially- 
covering circumnuclear gas, that both produces 
the time-invariant 'red wing' below the Fe K edge 
and modifies the spectrum above lOkeV. Such 
gas is already known in NGC 1365. Geometry- 
dependent Compton scattering must be taken into 
account, using 3D radiative transfer calculations, 
when making spectral models and when account- 
ing for the source total luminosity. The expected 
red wing, low Fe Ka equivalent width and high 
energy excess may be easily confused with 'reflec- 
tion' from an accretion disk, and spectral analyses 
to date have not established that there is any re- 
quirement for inner-disk reflection. 
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